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Abstract. It is known that the soft tail of the gamma-ray bursts’ spectra show excesses from the exact power-law 
dependence. In this article we show that this departure can be detected in the peak flux ratios of different BATSE 
DISCSC energy channels. This effect allows to estimate the redshift of the bright long gamma-ray bursts in the 
BATSE Catalog. A verification of these redshifts is obtained for the 8 CRB which have both BATSE DISCSC 
data and measured optical spectroscopic redshifts. There is good correlation between the measured and estimated 
redshifts, and the average error is Az « 0.33. The method is similar to the photometric redshift estimation of 
galaxies in the optical range, hence it can be called as ’’gamma photometric redshift estimation”. The estimated 
redshifts for the long bright gamma-ray bursts are up to z ~ 4. For the the faint long bursts - which should be 
up to 2 ~ 20 - the redshifts cannot be determined unambiguously with this method. 
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1. Introduction 

The ga mma-ray bursts (h e reafte r GRBs) of the long sub¬ 
group ( [Kouveliotou et ail , |1993[ ) detected by the BATSE 


instrument (Meegan et al 


2000| ) are at high reds hifts. The 


highe s t dire c tly measu r ed re dshift is at 2 = 4.5 ( Andersen 
et at, 200C|; Meszaros, 2001), but there are indirect con- 


isten c e of redshifts up to z 20 ([Meszaros fc Meszaro: 


1995, 1996; Horvath et al, 1996; Balazs et al, 1998). This 


ments ( |Klose| , P000| ) made possible to detect the afterglows 
and then the measurement of redshifts using optical spec¬ 


troscopy. The Current BATSE Catalog (Meegan et al 


siderations - based on BATSE data - predicting the ex- 


result is based on distribution densities and deals with the 
CRB redshifts in statistical sense only. This means that 
one may obtain the fraction of CRBs being at a given red¬ 
shift interval (see, for example, Schmidt (2001)), but one 
cannot obtain the redshift of a given CRB event. 

Th fre are only a few cases, when the obser~vration with 
the Be 


^000 ) consists of more than 1200 long bursts, but for only 
9 of them have redshift measurement (8 have redshifts 
and there is one GRB with an upper redshift limit). The 
unfortunate premature termination of the CGRO satellite 
prevents to increase this number further. There are other 
instruments observing « 100 bursts/year, but the typi¬ 
cal number of burst’s redshifts is only about a dozen/year 
(Meszaros, 2001). 

Hence, any method that could estimate the redshifts 
from X-ray/gamma-ray observations alone would be a 
great help. 

In Ramirez-Ruiz fc Fenimor^ (2000|) and Reichart 


ipoS AX s a tel l ite (|Piro et o/ | , o r o th er iu s tru- a linear relation between the intrinsic peak- 

- ^^ liiminnsitips nf CRBs and tbpir sn callpd ”variabilitips” 
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luminosities of GRBs and their so called ’’variabilities” 
was found. Similarly, Norris et al. (pOOOa ) found a relation 
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between the so called spectral lag and the peak-luminosity 
allowing to estimate the redshifts of long GRBs. These 
relations were calibrated on a few cases of GRBs, when 
GRBs were observed both by BATSE and other instru¬ 
ments measuring the optical redshift from afterglows. 
Then, having either the variabilities or the spectral lag of 
a given GRB, one can estimate its redshift. The physical 
meaning of the correlation between the variability (spec¬ 
tral lag) and the peak-luminosity remains unclear. These 
two methods can be combined (Schaefer et at, ^001) to 


determine redshifts if all the needed input parameters are 
available for the GRBs. 

In this article we present a new method of the estima¬ 
tion of the redshifts for the long GRBs. The situation is in 
some sense similar to the optical observations of galaxies, 
where the number of objects with broad band photomet¬ 
ric observations is much larger than the number of objects 
with measured spectroscopic redshifts. For galaxies and 


tion (Koo, 1985; Gonnolly et 

'll., 1995a; Gwyn & Hartwick 

, 1996 Sawicki et al. 

1997 

WangBahcall et al., 1998; 

Fernandez-Soto et al\. 

1999; 

Benitez, 200C; Gsabai et al., 

200 c; Budavari et al., 

200 c, 

2001) achieved a great sue- 


energy channels, the observed luminosity (in units pho¬ 
tons/s) for a burst at redshift z would be the following: 


|■{l + z)E2 /■(l + ZjJii 

/ L{E)dE, L 4.3 = / L{E)dE, (1) 


L2,1 — 




(1-\-z)E4 
{1 + z)E3 


For the observed fluxes (P) similar equations can be 
used. Let us define the following flux ratio: 


R{E 4 , E 3 , E 2 , El, z) — 


.^4,3 — ^2,1 -^4,3 — ^ 2,1 

La,3 + L 2 .I P4,3 + ^2,1 


( 2 ) 


which in general depends on the redshift, since L 2 ,i and 
L 4_3 depends on the redshift. 

Assume for the moment that one observes a pure 
power-law spectrum. This means that Le oa E~°‘ holds, 
where the exponent a is a real number. In this special case 
one could prove easily that for any redshift z 


rz=0 T z=0 

^ 4,3 ^ 2,1 

T Z—0 I T z—0 

-^4,3 -^2,1 

where 


La,3 — L2 
La ,3 + L 2 


cess in estimating redshifts from photometry only. Here we 
present a method that is quite similar to these methods; 
hence we call it as gamma photometric redshift estimation 
(GPZ for short). We utilize the fact that broadband fluxes 
change systematically, as characteristic spectral features 
redshift into, or out of the observational bands. Hence, 
contrary to the variability and spectral lag methods, this 
technique has a well defined physical meaning. 

The article is structured as follows. First, using a spec¬ 
tral model for GRBs we deduce an expected relation be¬ 
tween a measurable quantity (peak flux ratio) and the 
redshift (Sect. 2). Having this relation, we verify it on the 
existing sample of a few GRBs having measured redshifts 
(Sect. 3). Because both Sect. 2 and Sect. 3 suggest that 
this method is usable. Sect. 4 presents the estimated red¬ 
shifts for hundreds of long GRBs. In Sect. 5 we discuss 
and summarize the results. 


2. Gamma photometric redshift estimation 

To understand our method in this Section we outline the 
general scheme of broadband observations. The method 
is generally the same both for the optical and gamma- 
ray ranges. The only major difference is that in the X-ray 
and gamma-ray range the extra- and intergalactic medium 
have negligible effects, but the optical photons are atten¬ 
uated. 


L^=o= f \{E)dE, Llf^ = [ \{E)dE. 

J E^ £^3 


( 3 ) 


( 4 ) 


This means that in this special case R = 
R(Ea, E 3 , E 2 , El, z) is not depending on z. 

Of course, in the real situation, the spectrum has got 
a more complicated form, and hence R is depending on z; 
this will be the effect that we will use for redshift estima¬ 
tion. 

In addition, in the real situation, the incident spec¬ 
trum measured by the detector is convolved with the 
detector’s response fun ction defined by response matrix 
( Pendleton et al. , 1994 ) resulting the measured flux of the 
corresponding channel. For the channel with energy range 
E 2 > E > El the measured flux Pi ,2 is therefore given by 


Pi ,2 = [ ' P{E)c{E) dE, 

JE^ 


( 5 ) 


where c{E) is the detector’s response function. The similar 
holds for the second channel, too, with the same c{E). 
Hence, in general, R is depending both on the spectrum 
and the response function. 

If the rest-frame spectrum for a GRB is known, one is 
able to calculate the theoretical P as a function on z. Then 
these values can be compared with the flux ratio obtained 
from the broadband measurements {Rmeas)- The redshift, 
where (P — Rmeas)^ is minimal, could give the estimated 
gamma photometric redshift. 

Regarding this gamma photometric redshift estima¬ 
tion, the major problem comes from the fact that the 
spectra are changing quite rapidly with time; the typical 


Let us take two different instrumental channels defined timescale for the time variation is ~ (0.5 — 2.5) s (Ryde & 


by Ea > E 3 and E 2 > Pi. If one would know the rest- Svensson, 1999, 2000). Hence, if possible, one should con- 


frame energy spectrum (L{E)) of the burst, for a per¬ 
fect instrument that captures all the photons in the above 


sider spectra which are defined for time intervals smaller 
than this characteristic time. Therefore, we will consider 
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the spectra in the 320 ms time interval (i.e. in five 64 ms 
time intervals), with the peak-fiux being at the center of 
this time interval. 

In the following we will assume that the spectrum has 
the same shape around the time of the peak-fiux for all 
long bursts. Unfortunately we do not have any deep the¬ 
oretical or observational evidence for this assumption, in¬ 
stead we will test our assumption on GRBs, where spec¬ 
troscopic redshifts are available (next Section). Because, 
this assumption seems to be acceptable, in Sect. 4 we will 
use R to estimate z for long GRBs. 

3. Application on GRBs: Calibration 


a=3.2 p=0.5 £,,,=90 keV 



It is well-known (Band et al., 1994; Amati et ai, 2002) 


that the time-integrated average spectra of GRBs can be 
approximated by a broken power-law; the break is at some 
energy Eq. The typical rest-frame energy for Eo is above 
~ 500 keV (Preece et al, pOOOa ; Preece et al, 2000h| ), but 
this might vary for different GRBs. 

Of course this broken power-law spectrum is simply 
an approximation: first, because the break around Eg may 


have a more complicated form (Preece et at, 2000a; Preece 


et al\ , ^OOOb ), and, second, because at low rest-frame ener¬ 


gies (around ~ 80 keV) there may be essential departures 
from the power-law. This is the so called soft-excess, which 
is confirmed for ~ 15% of GRBs on the high confidence 
level (Preece et al, 1996, 2000a; Preece et al., 2000h| ); and 
for the remaining GRBs the soft-excess seems to occur, too 


( Preece et al , 1996| ). 


Based on this, we construct our template spectrum that 
will be used in the GPZ process in the following manner: 


Let the spectrum be a sum of the Band’s function (Band 


et al., 1994), and of a low energy power-law function tak¬ 


ing the form 

iovE<Eo L{E) =a{E/Egr)~°'+ a{E/Egr)~^ ( 6 ) 
for E>Eo L{E) = a^[E/Eg)-\ (7) 

where 03 = a[{Eo/Ecr)~°‘ + {Eo/Egr)~^] comes from the 
normalization. In this spectrum there are six parameters, 
but the amplitude a is for R unimportant, and need not 
be specified. We fix all the above parameters according 
to the available literature, so there are no free tunable 
parameters in our method. The low energy cross-over is 
at Ecr = 90 keV, Eg = 500 keV, and the spectral indices 
are a = 3.2, /3 = 0.5 and 7 = 3.0 (Preece et al 
20 ^. 


Preece et al. 


2000 a; 


Fig. 1. The detector response function c{E) and the be¬ 
haviour of a template spectrum at different redshifts. 


public BATSE Gatalog ( Meegan et at , 2000 ). During this 
time inter val the change of our tem plate spectrum is still 
negligible ( Ryde fc Svensson , 20001 ) . 

We have also checked the robustness of the PER 
against the integration time around the peak. Both the 
doubling of the integration time (for 640 ms) and its skew¬ 
ing around the peak did not change significantly the val¬ 
ues of PER. All this means that the template spectrum 
defined by Eqs.(6-7) seems to be a good approximation 
for the 320 ms time interval around the peak. This is in 


fact expectable from earlier studies of spectra (Band et al. 


1994; Ryde fc Svensson, 1999, 2000). 


The 4 energies for the BATSE instrument are: Ei = 25 
keV, E 2 = E 3 = 55 keV, E 4 = 100 keV. Using the de¬ 


tector response matrices (Pendleton et al, 1994) one can 


calculate the observed counts and flux for any incoming 
spectrum. In Fig. we show a typical response function 
c{E). The response function is different for each burst, but 
using the BATSE DRM data one can use the actual re¬ 
sponse function for every burst. Fig. also demonstrates 
the behaviour of the spectrum at different redshifts. Going 
from z = 0 to higher redshifts one can see that the soft- 
excess moves from the second channel to the first one and 
then leaves the range of this detector around z « 4. 

Before starting the detailed investigation of the fluxes 
that one can get using the template spectra and the re¬ 
sponse matrices by Eq.( 6 ), let us test the correctness of 
the template spectrum in a simple way. Let us introduce 
the peak flux ratio (PFR hereafter) in the following way: 


As we remarked above, the spectrum is rapidly chang¬ 
ing. But our assumption is that the spectrum has a char¬ 
acteristic shape around the instant of the peak. Now we 
have to chose a short time interval around the peak (max¬ 
imum of the total counts), during which the change of 
spectrum is still negligible, but the number of photons al¬ 
lows good signal-to-noise ratio. To be able to cut out such 
time interval around the peak-fiux, we need data with a 


PFR = 


hi — I 


12 


hi + h 


( 8 ) 


where fj is the BATSE DISGSG flux in energy channel 
Ei < E < Ej integrated for 320 ms around the peak 
flux. (I.e. five 64 ms intervals are summed - the middle 
one is where the flux is the biggest. Even during this time 


reason; ibly good time resolution. In our study we will use 


interval the change of spectrum is still negligible (Ryde 
& Svensson, 2000D). Our theory says that for the above 


the 64 ms resolution BATSE LAD DISGSG data from the template spectrum this ratio should increase with z. On 
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0=3.2 |3=0.5 E<,r=90 keV a=3.2 p=0.5 Ecr=90 keV 




Fig. 2. The theoretical PFR curves calculated from the 
template spectrum using the average detector response 
matrix. 


Table 1. The redshifts of GRBs that have both BATSE 
triggers and m e asure d spectro scopic redshif t s. Da ta com¬ 
piled by Klose ( 2000 ); see also Bloom et al. ( 2001 ). 


Burst 

BATSE 

trigger 

Z 

Remark 

980425 

6707 

0.00857 

SN1998bw 

970508 

6225 

0.8356 


970828 

6350 

0.9578 

no DISCSC data 

980703 

6891 

0.9676 


991216 

7906 

1.020 


990123 

7343 

1.6006 


990510 

7560 

1.6196 


971214 

6533 

3.4127 


980329 

6665 

3.5 

upper limit only 


Fig. 3. Experimental (points with error bars) and theoret¬ 
ical (dashed lines) peak flux ratio values for the 8 bursts 
that have both BATSE DISCSC data and measured red¬ 
shifts. 


0=3.2 p=0.5 Ecr=90 keV 


■ 


■ 


5N1998bw 

i i-"" 1, 

■ 4 

* —t 


upper limit 


■ 


0 1 2 3 4 5 6 


^measured 


Fig. we plot the theoretical PFR curves calculated from 
the above defined template spectrum using the average 
detector response matrices for the 9 bursts that have both 
BATSE data and measured redshifts. These bursts’ data 
are collected in Table 1. 

If we redshift the template spectrum and apply the 
corresponding detector response matrix of the given burst, 
we can get a PFR value for any redshift. Fig. shows the 
theoretical curves together with the seven PFR values cal¬ 
culated from observed GRB data (the DISCSC data for 
GRB 970828 are missing). We used the template spectrum 
defined by Eqs.(6-7) and the different response matrices 
corresponding to the observational conditions of the given 
burst. There is a clear trend: as expected from the above 
considerations, PFR increases with increasing redshifts up 
to 2 ~ 4. Except for the GRB associated with the super¬ 
nova, and for the GRB having only the upper limit, the 
remaining 6 GRBs have a clearly increasing PFR with 
increasing z. 

In the used range of z (i.e. for z^4:) the relation be¬ 
tween z and PFR is invertable. Hence we can use it to 
estimate the gamma photometrie redshift (GPZ) from a 


Fig. 4. The measured spectroscopic redshift values com¬ 
pared with gamma photometric redshift estimation for the 
8 considered GRBs. 


measured PFR. In Fig. ^the measured spectroscopic red¬ 
shifts are compared with GPZ values for 8 considered 
GRBs. The errorbars show the effect of counts’ Poisson 
noise only. 

Leaving out GRB associated with the supernova and 
GRB having upper redshift limit only, the estimation error 

is Az = « 0.33. 

In order to test the reality of the correlation between 
the soft excess and the redshift we made the null hypoth¬ 
esis that there is no relationship between these quantities, 
i.e the computed correlation is purely random. Assuming 
no true correlation between the soft excess and redshift 
the probability density of the computed quantity can be 
given by 


fix) 


1 r((jV-l)/2) 2 UN- 4)/2 

0Fr((iV-2)/2)^ ^ 
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Table 2. The linear correlation coefficients between 2 ®^®° 
and 


N 

Xc p 

= 1-/3 

Remark 

7 

0.9364 

0.9991 


6 

0.9655 

0.9991 

GRB 980425 excluded 

6 

0.7809 

0.9666 

GRB 971214 excluded 

5 

0.6618 

0.8881 

GRB 980425 -L GRB 971214 excluded 


where N is number of data points (Spiegel fc Stephens 


199E). Let Tc be the calculated correlation. Then the f3 


error probability at rejecting the null hypothesis is given 

by 


rc 

p =1- J f{x)dx. 

-r 

In Table 2. the N, Vc and the calculated level of signif¬ 
icance are shown for various cases. 

Although it seems that for GRB 971214 (where 

z = 3.4127) fits very well the estimation error without 
it is better: Az « 0.29. However the linear correlation 
coefficient here with N = 5 yields a much poorer Xc = 0.66 
with a p = 0.89 significance. 

We see that PFR (if calculable from observations for 
the given burst) is a quantity that may allow to determine 
redshift. Problems may arise from the fact that for any 
value of PFR two redshifts are possible - either below or 
above z ~ 4 (see Fig. ||), further measurements are needed 
to exclude one of the redshift. 


4. Application on GRBs: Estimation of the 
redshifts 


To avoid the problems with the instrumental threshold 
we exclude the faintest GRBs from the BATSE data. 


Simila rly to Pendleton et al. ( 1997 ) and Balazs et al 
( |l998| ) these events have a F 256 peak-flux (i.e. on 256 ms 
trigger scale) smaller than 0.65 photon/(cm^s). These 
GRBs are not discussed in this article. 

Further restriction comes from the fact that short 


GRB s are today taken a s differe nt phenomena ( Horvath 
et al , 200 c; Norris et al , 2000^. In addition, due to in¬ 


strumental effects ( Piro et al. , 2002| ), no spectroscopic red¬ 
shifts are known for this subgroup of GRBs. Hence, we do 
not apply our method for short GRBs. 

The reality of the intermediate subgroup of GRBs 
(Horvath, 199^; Mukherjee et al., 1998; Hakkila, et al 


200C; Balastegui et al., 2001; Rajaniemi et al, 2002; 
Horvath, 2002|) havin g remarkable sky angular distribu¬ 
tion (Meszaros et al, 2000a,Litvin et al, 2001) is un¬ 
clear yet. In any case, no spectroscopic redshifts are known 
also here. Hence, we exclude this subgroup, too. 


The refore, we restr i ct ou rselves to long GRBs defined 
by Tqo ( [Meegan et al , 2000| ) with T 90 > 10 s . There are 
1241 GRBs in BATSE Catalog fulfilling this condition. 
Deleting GRBs having no F 256 and having ^256 < 0.65 
photon/(cm^s) 838 GRBs remain. This sample is studied 
here. 


Introducing an another cut F 256 > 2.00 photon/(cm^s) 
we can investigate roughly the brighter half of this sam¬ 
ple. We will discuss the sample F 256 > 2.00 photon/(cm^s) 
(’’bright half’ sample having 343 GRBs) and F 256 > 0.65 
photon/(cm^s) (’’all” sample having 838 GRBs), respec¬ 
tively. 

As the soft-excess range redshifts out from the BATSE 
DISCSC energy channels around z « 4, the theoretical 
curves converge to a constant value. For higher z it starts 
to decrease. This is where the power-law breakpoint {Eg) 
is redshifts into soft energy range. This means that the 
method is ambiguous: for the given value of PFR one may 
have two redshifts - below and above z « 4. Because for 
the bright GRBs the values above z « 4 are practically 
excluded, for them the method is usable. In other words, 
using only the 25—55 keV and 55—100 keV BATSE energy 
channels, this method can be used to estimate GPZ only 
in the redshift range z P A; outside of this region the z 
vs. PFR relation is non-invertable (see Fig. ||). For high 
redshifts (above z « 4) the method gives two possible 
values. For faint GRBs the estimation also usable (at least 
in principle), but one has to decide by other arguments 
that either the redshift below z « 4 or above z « 4 is the 
correct value. 

Let us assume for a moment that all observed long 
bursts, we have selected above, have z < 4. Then we can 
simply calculate the z^^^ redshift for any GRB, which has 
calculable PFR from BATSE DISGSG data. Fig. || shows 
the distribution of the measured PFRs of the long GRBs 
having DISCSC data. The fact that the number of objects 
beyond the minimal and maximal theoretical PFR values 
(« —0.15 and « 0.37, respectively) is relatively small, is 
reassuring. 

Fig. H shows the distribution of the estimated derived 
redshifts under the assumption that all GRBs are below 
z « 4. The distribution has a clear peak value around 
PFR « 0.2, which corresponds to z « (1.5 — 2.0). 


5. Discussion 


Having the estimated redshifts shown on Fig. ^ one may 
ask: Are these redshifts really correct? 

There can be two different problems here. First of all, 
the method is based on the assumption that around the 
peak flux, the spectrum is the same for all the selected 
long GRBs. Second, the method gives degenerate result, 
with two possible redshift values. 

Concerning the first problem we could just hope that 
in the near future some theoretical or experimental ev¬ 
idence will confirm our assumption, but the situation is 


not worse than in the studies of Norris et al. (2000a) and 
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Fig. 5. The PFR distribution of the long GRBs having 
DISCSC data. 



Fig. 6. The distribution of the Gamma Photometric 
Redshift estimators of the long GRBs having DISGSG 
data. The distribution of the bright half of the BATSE 
Catalog is also shown. 


Reichart et al. ( pOOl ). These articles also suggest that de¬ 
spite the deeper understanding of the underlying physics, 
the procedure itself is usable. In addition, here the PFR-z 
relation is well supported by earlier independent observa¬ 
tions. 

Concerning the second problem we think that the great 
majority of values of z obtained for the bright half are cor¬ 
rect. This opinion may be supported by three independent 
arguments. First, the obtained distribution of GRBs in z 
for the bright half on Fig. | is v ery similar to the obtained 
distribution of ^chmidt (2001) (see Fig. 6 of that arti¬ 


cle). The luminosity-based redshift distribution (Schaefer 
et aZ.JlOOl) also suggest an uniformly rising GRB density 
out to z « 5. Second, as z moves into z^4 regime for the 
bright GRB, one would obtain extremely high luminosi¬ 
ties. Using Eq.l2 of Meszaros & Meszaroj (1996), there is 


a lower limit for the isotropic luminosity of the GRBs a 
value ~ lO^^ergs/s. (Note here that the precise value is, of 
course, calculable and is depending on the chosen cosmol¬ 


Fig. 7. The redshift distribution of the 17 GRBs’ with 
known redshift and the distributions from the Gamma 
Photometric Redshift estimators for different peak flux 
groups. 


ogy model and on the typical energy of emitted photons. 
For the purpose of this article this approximate value is 
enough.) This is an unacceptable high value for a lower 
limit, because typical luminosities a r e 1 0 ^^~^^erffs/s. 
(see, e.g.. Table 1 of Reichart et al. ( 2001 )). We cannot 
exclude that a few cases from bright GRB are at z^4, 
but - we think - in the bright half this cases are rare. 
Thirds, as an additional statistical test we compared the 
redshift distribution of the 17 GRB with observed red¬ 
shift with our reconstructed GRB z distributions (limited 
to the z < 4 range). For the F 256 > 0.65 photon/(cm^s) 
group the Kolmogorov-Smirnov test suggests a 38% prob¬ 
ability, i.e. the observed N{< z) probability distribution 
agrees quite well with the GPZ reconstructed function. 
Although the observed distribution suffers from strong se¬ 
lection effects this fact is nevertheless reassuring. 

For the faint GRBs being between F 256 = 0.65 
photon/(cm^s) and F 256 = 2.00 photon/(cm^s) the sit¬ 
uation is different. From Fig. ^ it follows that for z < 1.7 
GRBs should be dominated by faint objects. From this 
Figure one would obtain that GRBs are in average at 
smaller redshifts. This is clearly a wrong conclusion, which 
is caused by the false assumption that z < 4 for all the 
faint GRBs. We think that the majority of faint GRBs 
z should be changed into the value z^4. Unfortunately, 
we are not able to say, concretely which GRB has a great 
(zj^4), and which GRB has still a small (z^A) redshift. In 
addition, for these faint GRBs also the error of estimated 
redshifts should probably be bigger than Az « 0.33. 
Simply, we conclude that in the current form with the 
current data, our method is not applicable for the faint 
GRBs. 

The results of this work may be summarized as follows. 


1. Based on earlier observations of GRB spectra it is 
shown that the peak flux ratio (PFR) should be a well 
defined function of z. 
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2. The estimated redshifts from PFR are in good accor¬ 
dance with the known redshifts of the few GRBs in 
BATSE Catalog having spectroscopic redshifts 

3. All this allows us to calculate the redshifts of long 
GRBs. Unfortunately, due to the twofold character of 
the PFR curve, the method is usable only for bright 
GRBs. 

4. Redshift distribution of 343 bright long GRBs are de¬ 
termined (Fig. H). 
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